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ABSTRACT 

Context. Single-dish observations of molecular tracers have suggested that both star formation and an AGN can drive the gas chemistry 
of the central ~kpc of active galaxies. The irradiation by UV photons from an starburst or by X-rays from an AGN is expected to 
produce different signatures in molecular chemistry, which existing data on low-/ lines cannot distinguish, as they do not trace gas at 
high temperature and density. Depending on the angular scale of a galaxy, the observed low-/ lines can be dominated by the emission 
coming from the starburst ring rather than from the central region. 

Aims. With the incorporation of high-/ molecular lines, we aim to constrain the physical conditions of the dense gas in the central 
region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. 
Methods. We used the James Clerk Maxwell Telescope to observe the /=4-3 transition of HCN, HNC, and HCO+, as well as the CN 
Nj = 25/2 - I3/2 and Nj = 35/2 - 25/2, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line 
intensity ratios and radiative transfer models. We discuss the results in the context of models of irradiation of the molecular gas by 
UV light and X-rays. 

Results. A first-order estimate leads to starburst contribution factors of 0.58 and 0.56 for the CN and HCN 7=1-0 lines, respectively. 
We find that the bulk emission of HCN, HNC, CN, and the high-/ HCO^ emerge from dense gas (n(H2) > 10^ cm"^). However, the 
low-/ HCO^ lines (dominating the HCO^ column density) trace less dense (n(H2) < 10^ cm"^) and colder (Tk < 20 K) gas, whereas 
the high-/ HCO+ emerges from warmer (> 30 K) gas than the other molecules. We also find that the HNC/HCN and CN/HCN line 
intensity ratios decrease with increasing rotational quantum number /. 

Conclusions. The HCO^ /=4-3 line intensity, compared with the lower transition lines and with the HCN /=4— 3 line, support the 
influence of a local XDR environment. The estimated A'(CN)/A'(HCN) ~ 1 - 4 column density ratios are indicative of an XDR/AGN 
environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks. 

Key words, galaxies: ISM — galaxies: individual: NGC 1068 — galaxies: Seyfert — radio lines: galaxies — radio lines: ISM — 
ISM: molecules: HCN, HNC, CN, HCO+ 



' 1. Introduction The prototypical Seyfert 2 galaxy NGC 1068 (located at 

a distance ~ 14.4 Mpc and with a bolometric luminosity of 
^ Active galaxies often have high concentrations of molecular gas ^ lO^'L©) is a known source where the star formation and/or 
in their central (<kpc) region, which may be feeding a central activity can be debated as the main driving mechanism 
active galactic nucleus (AGN) or an starburst (e.g. Sanders & ^f its gas chemistry. It hosts a circumnuclear starburst ring of 
Mirabel 1996; Gao & Solomon 1997, 1999). Molecular gas ir- ^^out 2.5 - 3 kpc in diameter (Schinnerer et al. 2000), which 
radiation by hard X-ray photons (> 1 keV) emitted during the delimits a ~ 2.3 kpc stellar bar (ScoviUe et al. 1988). The de- 
accretion process in the proximity of AGNs, leads to an X-ray ,tro„g emission coming from the starburst ring 
dominated region (XDR) (e.g. Maloney et al. 1996; Lepp and indicates that its massive star formation activity is being fed by 
Dalgarno 1996; Meijennk & Spaans 2005). In an starburst re- ^^^^ amounts of molecular gas (e.g. Planesas et al. 1989, 1991; 
gion, the formation/destruction of molecules is instead driven by ^gjfgj & Blitz 1995; Schinnerer et al. 2000). However, signifi- 
far-UV photons (6 - 13.6 eV) emitted by O and B stars, which ^^^^ emission also arises from a circumnuclear disk (CND) 
leave the fingerprint of a photon-dominated region (FDR) (c.f. ^^out 200 pc in diameter (Solomon & Barret 1991; Schinerer 
Hollenbach & Tilens 1999, and references therein). An increased 200O). Earlier observations of high HCN/CO luminosity 
X-ray ionization of molecular clouds can enhance the abundance ^^^^^^^ measured in the CND of NGC 1068, suggested an abnor- 
of several molecules (e.g. CN, HNC, HCN and HCO+) with re- jjial (by Galactic standards) chemistry of molecular gas (Tacconi 
spect to their abundances found in PDRs (e.g. Lepp & Dalgarno ^9941 Sternberg et al. imi Later observations of HCN, 
1996; Meijennk & Spaans 2005; Meijennk, Spaans & Israel ^^0+ and HCO, among other molecules (Usero et al. 2004}, 
2007). Conversely, the chemical state of the molecular gas is a ^^ggg^j j^at the chemistry in the CND of NGC 1068 cone- 
tracer of the dominant energy source of AGNs. sponds to that of an XDR. On the other hand, single-dish obser- 
vations of other high-density tracers (HNC and CN) in the CND 
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of NGC 1068 show that the HNC/HCN and CN/HCN line ra- 
tios favor a PDR, rather than an XDR scenario (Perez-Beaupuits 
et al. 2007 ). This independent observational evidence for both 
PDR and XDR scenarios is consistent with a model that consid- 
ers an AGN and starburst component to be present in NGC 1068 
(Spinoglio et al. l2005l l. However, a later estimate indicates that 
the age of the most recent episode of star formation, within the 
central 10-100 pc, is about 250 + 50 Myr (Davies et al. 2007). 
This result suggests a recent, although no longer active, starburst 
activity in the CND of NGC 1068. 



Nevertheless, it was noted in Perez-Beaupuits et al. ( I2007I I 
that the result of a favored PDR scenario in the CND of 
NGC 1068 may be misleading, since the AGN contribution 
(through XDR eflFects) is typical of an smaller angular scale 
than the starburst contribution (PDR). X-rays are known to pen- 
etrate deeper than UV photons into a single molecular cloud. 
But its effects are rather local, since the energy flux of the ra- 
diation decreases with radial distance from the source. On the 
other hand, an starburst region is more spread, and several O/B 
stars will typically be closer to a molecular cloud than a single 
X-ray source. Thus the PDR effects will dominate the bulk of 
the gas. Hence, the smaller angular size of the action zone of 
an X-ray source (the region around an X-ray source that domi- 
nates over UV driven chemical effects) will be affected by beam 
dilution, and the radiation emerging from the starburst compo- 
nent can dominate the overall emission collected by a single dish 
beam. This is a known drawback of single dish observations of 
extragalactic sources that, depending on the diameter of the tele- 
scope used, complicates the efforts to disentangle the AGN ac- 
tivity from that of the starburst ring, when considering the in- 
tensities of the lower transitions of molecular lines. Given that 
the starburst contribution is likely dominant at radii R ~ 14" 
in NGC 1068, but not in the CND at R ~ 2", the low-7 lines 
observed with single dish telescopes (with beam > 25") can be 
dominated by the contribution from the starburst ring, and if un- 
corrected, can be misleading. 



In this paper we present observations of high-7 lines of HCN, 
HNC, CN and HCO+ on NGC 1068. Due to the smaller beam 
size, these lines are less affected by the contribution from the 
starburst ring, than the 7=1-0 lines. We also present a first or- 
der estimate of the starburst contribution, based on independent 
observations of the 7=1-0 lines, and correct the intensities of 
these lines accordingly. Since HNC is the isomer of HCN, the 
radiation of these molecules is expected to arise from gas with 
the same physical conditions. However, radiative transfer models 
based on the yr^ line ratios of these molecules showed that the 
HNC emissions arise from less dense gas than HCN, when con- 
sidering a kinetic temperature of about 80 K (Perez-Beaupuits 
et al. I2007I I. We compare these previous results with those ob- 
tained from the new line ratios. The molecular abundances 
estimated from the radiative transfer models are subsequently 
compared to those obtained with our PDR and XDR models. 



In §2 we describe the observations. The results (spectral lines 
and line intensities) are presented in §3. The modelling and anal- 
ysis of the excitation conditions and abundances of HCN, HNC, 
HCO^ and CN are presented in §4. The conclusions and final 
remarks of this work are presented in §5. 



Table 1. Beam sizes & efficiencies. 



Molecule 


Transition 


V 


HPBW 






J/Nj 


[GHz] 


["] 




CN 


25/2 - I3/2 


226.875 


21 


0.67 


CN 


35/2 - 25/2 


339.517 


15 


0.63 


HCN 


4-3 


354.505 


14 


0.63 


HCO+ 


4-3 


356.734 


14 


0.63 


HNC 


4-3 


362.630 


14 


0.63 



2. Observations and data reduction 

We have used the 16-pixel receiver HARP (Heterodyne Array 
Receiver Programme) on JCMlQ, during December 2006, to ob- 
serve the 7=4-3 Hnes of HNC (363 GHz), HCN (355 GHz) 
and HCO+ (357 GHz), as well as the CN Nj = 35/2 - 25/2 
(339 GHz), toward the center of the Seyfert galaxy NGC 1068 
(ra.= 02H2'"40.7^ dec.= -00°00'48"). Beamsizes and efficien- 
cies are shown in Table [1] We used HARP as a single pixel in 
staring mode, since it was the only receiver available to cover 
the high frequencies. The 10* pixel was used to point the tele- 
scope at the nuclear region of NGC 1068. Detection in the other 
pixels was negligible and beyond the scope of this work, since 
they point outside the circumnuclear disk. Pointing was checked 
regularly on SiO masers and the RMS was found to be lower 
than 2". Since HARP was not fully commissioned at the time 
of the observations, we regularly checked the calibration of the 
receiver observing the CO 7=3-2 line in the reference source 
CRL168 and we compared it with the representative spectrum 
obtained with the former receiver B3. It was found that the cal- 
ibration CO 7=3-2 line observed with HARP was fainter than 
the reference spectrum, and our data was corrected according 
to the variations of the calibration observed during the run. We 
also used the receiver A3 to observe the CN Nj - 2^/2 - I3/2 
(226 GHz) line. The backend was the ACSIS correlator, and 
when observing with HARP we used a bandwidth of 1 GHz 
since the 1.8 GHz bandwidth was showing some anomalies in 
the overlap area and a higher system temperature. Unfortunately, 
this limited our baseline coverage. Nevertheless, we used 2048 
spectral channels which gives a resolution of about 0.4 km s"' in 
the original spectra. We obtained different numbers of scans (on- 
source observations) for all the observed lines, with integration 
times between 5 and 10 minutes each, to allow for calibration 
observations between the scans. 

After a weighted (by system temperature and integration 
time) average of the scans, a baseline was removed from the av- 
erage spectra and the resolution was lowered to reduce the level 
of noise. A first-order polynomial was used in most of the spectra 
to fit the baselines. Only for CN Nj - 35/2 - 25/2 a second order 
polynomial was needed. The velocity resolution was reduced to 
20 kms which represents less than 10% of the line widths. 
Small variations in the channels chosen to define the baselines 
in the original averaged spectra, lead to variations of < 6% in 
the integrated intensities of the final spectra shown in the next 
section. The average system temperature was about 300 K for 
total on-source integration times between 35 minutes and 4.1 
hours. The S/N levels (the peak antenna temperatures, obtained 
from the gaussian fit, with respect to the rm.s. of the baselines) 



' The JCMT is operated by the Joint Astronomy Center on behalf of 
the Science and Technology Facilities Council of the United Kingdom, 
the Netherlands Organization for Scientific Research, and the National 
Research Council of Canada. 
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Table 2. NGC 1068 line parameters. 



1000 1300 
Velocity [km/s] 



HCN /=4-3 




800 1000 1300 1400 

Velocity [km/s] 



; - HNC J"=4-3 




1000 1200 
Velocity [km/s] 




Velocity [km/s] 



Fig. 1. Molecular line emissions in NGC 1068. The velocity resolution 
was reduced to 20 km s~'for all the lines. The spectra are centered with 
respect to the heliocentric systemic velocity v^ys = 1 137 km s~' . 



Molecule-y 


V 

[km s"'] 


n 

[mK] 


[kms"'] 


CNA'y = 25/2-13/2 


1151 + 15 
1444+42 


21.4+ 1.7 
7.1+ 2.1 


265+35 
265 " 


CNTVj = 3j/2-25/2 


1252+33 


2.2+ 0.8 


132" 


HCN J=4-3 


1123+ 6 


41.5+ 2.4 


202+14 


HNC y=4-3 


1116+18 


6.7+ 1.1 


218+47 


HCO+ 7=4^3 


1112+ 8 


14.9+ 1.4 


166+19 



The uncertainty of tliis parameter was larger tlian 100%, if let free. 
So we set its value to the one found for the main gaussian compo- 
nent. 



Table 3. Velocity-integrated intensities. 



Molecule 


Transition 


Intensity 




J/Nj 


[Kkms^'l 


CN 


25/2 - I3/2 


11.7+ 1.4 


CN 


35/2 - 25/2 


1.6+ 0.3 


HCN 


4-3 


13.9+ 1.6 


HNC 


4-3 


2.7+ 0.3 


HCO+ 


4-3 


3.8+ 0.5 



a) Total velocity integrated temperature considering all the compo- 
nents of each spectral line. The values correspond to the main-beam 
brightness intensity, I,„ij = j T'^lrj^^hSv, in [K km s"']. 



ranged between 3 and 7. The software package XS (written by 
P. Bergman, Onsala Space Observatory) was used to reduce the 
data and fit the Gaussians. 



3. Results 

The observed molecular line emissions are shown in Figure [T] 
and the line parameters of the gaussian fits are summarized 
in Table |2] The spectra are centered with respect to the he- 
liocentric systemic velocity Vsys - 1137 kms"', taken from 
the NASA/IPAC Extragalactic Database (NED). The velocity- 
integrated intensities are shown in Table [3] 

The HCN 7=4-3 spectrum shows the same single peak struc- 
ture as the lower-7 transitions reported by Perez-Beaupuits et 
al. (2007 ) and Curran et al. (f2000l l. The average line width is 
about 25% (~ 70kms"') narrower than the average width ob- 
served in the lower 7- transitions. 

The double peak structure seen in the HNC 7=1-0 and 7=3- 
2 spectral lines by Perez-Beaupuits et al. (2007 ) is not present in 
the 7=4-3 line. This may be due to the smaller beam size and 
to the lower signal to noise ratio of our data. The center velocity 
of the gaussian fit is just about 10 km s"' away from the center 
velocities of the main component in the lower transition lines 
(within ~ Icr). Since this is a single component line, its average 
width is about 50 km s ' wider than the line widths of the lower 
transition lines reported by Perez-Beaupuits et al. (I2007I I, but 
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similar to the line width of the HNC 7=1-0 line reported by 
Hiittemeister et al. ( |1995l l. 

The HCO^ spectrum shows a single peak structure, and ap- 
pears to be narrower than the 7=1-0 and 7=3-2 lines reported 
by Kjips et al. ( |2008| ). Our 7=4-3 line has an average line width 
of about 166 km/s, while the lines reported by Krips et al. have 
average line widths of about 240 km/s. This may be an indication 
that the 7=4-3 line traces a different gas than the lower transi- 
tions. But due to the relatively low S/N~ 6, deeper observations 
are needed. 

Using receiver A3, the two CN spingroups at 226.8746 GHz 
(AT = 2-1, 7 = 5/2 - 3/2, F = 7/2 - 5/2) and 226.6596 GHz {N 
= 2-1,7 = 3/2 - 1/2, F = 5/2 - 3/2) were detected. The central 
velocity of the second spingroup is about 17 kms ' away from 
the expected velocity, which is within the Icr error. In the opti- 
cally thin limit, the second spingroup is expected to be a factor 
< 0.56 weaker than the main spingroup. The factor we obtained 
from the gaussian fit in Figure [T]is ~ 0.33 + 0.1, which means 
that this line is close to optically thin. 

Observing CN with HARP at ~ 339.5 GHz, we also detected 
a double peak structure. However, the two spingroups at 339.517 
GHz (N = 3-2, 7 = 5/2- 5/2, F = 7/2 - 7/2) and 339.476 GHz 
(F = 5/2 - 5/2) are severely blended. The separation between the 
two peaks is 176 + 36 kms"', whereas the separation between 
the farthest spingroup (at 339.447 GHz, F = 3/2 - 3/2) that we 
could detect in our bandwidth is expected to be at just about 62 
km s"' towards the higher velocities. So we do not think the dou- 
ble peak structure is due to another spingroup. Nevertheless, two 
gaussian components are needed to fit the line profile. Given the 
noise in the spectrum, the line width of the second component 
was set to the fine width found for the main component. 

4. Discussion 

4.1. Structure of the CND and limitations of the model 

Millimeter/sub-millimeter and mid-IR high spatial resolution 
observations of the nuclear region of NGC 1068 support the 
classical picture of a rotating torus (e.g., Tacconi et al. 1994; 
Schinnerer et al. 2000; Galliano et al. 2005; Poncelet et al. 
2007). However, in some previous works (e.g., Schinnerer et al. 
2000) it was speculated that the classical (large-scale) molecu- 
lar torus is not always needed to explain the occultation of the 
nuclear regions. Recent mid-IR observations have shown the ex- 
istence of two molecular streamers feeding the galactic nucleus 
of NGC 1068 (Tomono et al. 2006; Sanchez et al. 2009). These 
observations suggest that the occultation of the nuclear region 
can also be associated with a nuclear concentration of molecu- 
lar gas and dust that forms an optically thick amorphous clumpy 
structure (Sanchez et al. 2009). Although it is very likely that 
this structure encloses smaller infalling clouds, they are not re- 
solved at the 0.075" (scale size ~ 5 pc at a distance of 14.4 Mpc) 
resolution of the Very Large Telescope (VLT). 

In the case of millimeter and sub-millimeter interferom- 
eter maps of high density tracers (e.g., Tacconi et al. 1994; 
Garcia-Burillo et al. 2008), not even the streamers are resolved. 
This represents a limitation of the current millimeter and sub- 
millimeter observations in comparison with visual and IR wave- 
length ranges, with which the morphology and dynamics of the 
gas can be studied in more detail. This limitation is even more 
severe in single dish observations for which, depending on the 
lines observed, the emitting regions (or clouds) can be of a 
much smaller scale (~ 1 order of magnitude) than the beam 
size of the telescope. In particular, our single dish observations 



of HCN, HNC, CN and HCO+ are sensitive to the high den- 
sity (~ 10^ cm"^) clouds that can be embedded even inside the 
feeding streamers observed by Sanchez et al. (2009). The size of 
these dense clumps is expected to be from < 1 pc up to a few 
pc, embedded in large molecular clouds of the order of tens of 
parsecs. In the near future, only the Atacama Large Millimeter 
Array (ALMA) is expected to resolve these smaller clouds at 
millimeter and sub-millimeter wavelengths with resolutions (at 
1mm) between 0.1" and 0.01" (scale size 7 - 0.7 pc at the dis- 
tance of NGC 1068). 

Ideally, we would like to be able to observe and to model 
the physical environment in extra-galactic sources, just like it 
has been done for the Milky Way. Relatively complex radiative 
transfer codes like, for instance, RATRAN (Hogerheijde & van 
der Tak 2000) and ^3D (Poelman & Spaans 2005), can be used 
to model the internal dynamics, temperature and density struc- 
ture of individual molecular clouds. However, the detailed struc- 
ture and signatures of the individual clouds are convolved and 
smeared out with a large (> 10") beam in single dish extra- 
galactic observations. Therefore, we do not have the informa- 
tion needed to constrain the parameters related to the detailed 
internal structure of the clouds. What we have is the convolved 
intensity and line profile of the emission emerging from all the 
clumps that fall within the telescope beam. Therefore, all model- 
ing efforts are limited to estimate only the bulk of the gas prop- 
erties. For this purpose, the less complex radiative transfer code 
RADE}<Q (Van der Tak et al. I2007I I is a sensitive model for es- 
timating these average physical conditions of the high density 
gas from single dish observations towards extra-galactic sources. 
However, due to the lack of spatial resolution in our data, and the 
limitations of our model, the following discussion and the con- 
clusions derived from them should be taken as qualitative orders 
of magnitude. Higher spatial resolution, such as the resolutions 
that ALMA will provide, are necessary to reach definitive con- 
clusions. 

4.2. Line ratios and starburst contribution 

The observed line ratios of molecular transitions have been 
largely used as a tool for estimating the physical conditions of 
the molecular gas. However, the conversion from to Tr is 
uncertain through the source size of the emitting region and, 
hence, the derived physical conditions {Tg, n{H2) and A^) depend 
strongly on this conversion. Special care regarding this issue has 
been taken in previous work (e.g. Tacconi et al. 1994, Usero et 
al. 2004, Perez-Beaupuits et al. 2007, Krips et al. 2008). 

The high resolution map of the CO 7=1-0 Hne in NGC 1068 
reported by Schinnerer et al. (2000) shows that the morphology 
of this galaxy consists mostly of a nuclear region (circumnu- 
clear disk) and spiral arms (starburst ring). Since our aim is to 
use the HCN, CN, HNC and HCO+ data as a diagnostic tool for 
AGN-like environments, we are interested only in the emission 
emerging from the CND of NGC 1068. However, the starburst 
ring could potentially be contributing to the observed emission 
and, hence, influence the results of the radiative transfer models, 
as noted by Perez-Beaupuits et al. (2007). This is particularly 
relevant for the 7=1-0 lines of all the molecules because of the 
larger beam size. Therefore, the line intensity ratios, involving 
the 7=1-0 transitions, could be misleading. 

A priori, the most straightforward way to estimate the con- 
tribution from the starburst ring would be to use interferome- 
ter data, for instance, from the HCN 7=1-0 map reported by 

" http://www.sron.rug.nl/~vdtak/radex/radex_manual.pdf 
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Tacconi et al. (1994). Although this is possible, using interfer- 
ometer data is somewhat delicate due to the missing flux. Even 
though the zero spacing problem can be resolved through spec- 
tra observed with a single-dish telescope in the direction of the 
interferometer phase center, only part of the short spacing infor- 
mation can be recovered. Besides, high quality interferometer 
and single-dish data are required in order to minimize the ar- 
tifacts in the final combined data. Because the uncertainties in 
this technique are higher, and because high quality data are not 
available for all the lines needed, we derive below a first-order 
estimate of the starburst contribution, based on single-dish data 
only. We use the available 7=1-0 lines of HCN, HNC and CN 
observed with diff'erent telescopes (SEST, OSO, IRAM 30m). 

The total flux or intensity Ftotal, measured with a beam 
of FWHM=6'mb that includes emission from the starburst ring 
(Fsb), of average radius r ~ 14", and from the CND (Fcnd, with 
an smaller radius of about 1"), can be decomposed as follows: 

■FtOTAL = -FCND * WCND + -FsB * WSB, (1) 

where wcnd and wsb are weighting functions that can be esti- 
mated as Wcnd ~ 1, because the size of the CND is small com- 
pared to the beams of any of the available single dish data, and 



Table 4. Line intensity ratios between the peak antenna temper- 
atures, assuming a source size of 1.5". 



Molecule 


Transitions 


Ratio " 


CN 


25/2-I3/2 

'3/2 ^^1^2 


0.61±0.39 


CN 


35/2-25/2 
25/2-I3/2 


0.13±0.04 


HCN 


3-2 
1-0 


1.01±0.53 


HCN 


4-3 
3-2 


0.49±0.14 


HNC 


3-2 
1-0 


0.38±0.18 


HNC 


4-3 

3-2 


0.30±0.10 


HCO+ 


3-2 
1-0 


0.12+0.04'' 


HCO+ 


4-3 
3-2 


1.37+0.46 



a) See Perez-Beaupuits et al. (2007) for the 7=1-0 and 7=3-2 tran- 
sitions. The low-7 lines of HCO^ were obtained from Rrips et 
al. (2008). 

b) This ratio slightly differs from that found by Krips et al. (2008) be- 
cause we use the ratio between the peak main beam temperatures 
and we corrected the HCO^ 7= 1-0 line for the starburst contribu- 
tion factor. 



Wsb ~ exp 



'14" 

-4 * ln{2) * I 

1 ^mb 



(2) 



In eq. (2) we assume that the contribution comes only from 
an axisymmetric ring-like distribution of radius r - 14", ne- 
glecting the contribution from the bar or any other feature at 
intermediate radii. Since we have at least two independent ob- 
servations of the 7=1-0 line, with diff'erent telescopes of beams 
9mh\ and 0nib2, we can have a first-order estimate of the starburst 
contribution as: 



WSBl - WSB2 



(3) 



Hence, from eq. (1), the starburst contribution factor /sb of 
a line observed with any telescope can be estimated as: 



F5B * Wsb 



SB 



TOTAL 



(4) 



So we can estimate the intensity that emerges from the CND 
from the observed total intensity as Fcnd ~ -Ftotal * (1 - /sb)- 
This means that, the larger /sb, the larger will be the |^ ratios. 
Note that, from equation (2), /sb depends on the beam size 6^\, 
of the telescope used for the observations. Using the peak in- 
tensities from Perez-Beaupuits et al. (2007) and from available 
IRAM 30m data for the same CN, HCN and HNC 7=1-0 lines 
(Usero et al. 2009) and the corresponding telescope beams, we 
estimate the /sb factors as 0.58+0.25 for CN (SEST, 9^^ = 45"), 
0.56 + 0.32 for HCN (OSO, fl^b = 44") and 0.99 + 0.50, for HNC 
(SEST, 0„b = 55"). 

The similar (within ~ 4%) starburst contribution factors 
found for CN and HCN, as observed with SEST and OSO, re- 
spectively, are as expected since the respective beam sizes are 
also similar (within ~ 2%). Instead, the high contribution fac- 
tor estimated for HNC would indicate that only 1 % of the ob- 
served emission emerges from the nuclear region, in contrast to 
the 42% and 44% found for CN and HCN. If this was real, we 
would expect to see a pronounced double peak structure, with 
the proper velocity centers, in the HNC 7= 1-0 line profile. Since 



this is not the case, we think that the substantial difference be- 
tween the starburst contribution factors of HNC and the other 
molecules might be due to some problems (e.g., pointing or cal- 
ibration errors) in the spectral line reported in Perez-Beaupuits 
et al. (2007), which was noted to be different than the one pre- 
viously reported by Hiittemeister et al. (1995). Thus, we use the 
HNC 7=1-0 line from Hiittemeister et al. for the following anal- 
ysis. However, this line, as well as the low-7 HCO^ lines re- 
ported by Krips et al. (2008), was observed with the IRAM 30m 
telescope only. Therefore, we cannot estimate the corresponding 
starburst contribution factor, due to the lack of independent ob- 
servations. Hence, for HNC and HCO^ 7=1-0 lines we adopt 
the value /sb = 0.45 + 0.31 estimated from the HCN observa- 
tions done with the 30m telescope (0n,b - 28"). The starburst 
contribution factors for the higher-7 transitions are expected to 
be negligible because of the smaller beams at higher frequencies, 
so we can assume these to be zero. 

After correcting the 7=1-0 lines for the starburst contri- 
bution, we correct the peak intensities for beam dilution as 

The ratios between the line intensities are 



T' 



shown in Table |4] We use the ratios between the peak inten- 
sities of the main components in the Gaussian fit shown in 
Table |2] The peak intensity corresponds to the emission emerg- 
ing from regions with velocities close to the systemic velocity of 
NGC 1068, and the ratios obtained from these values should be 
similar to those obtained from the velocity-integrated intensities, 
if the line widths of the corresponding components of the diff'er- 
ent transitions are also similar. The peak intensities of the low-7 
transitions of HCN and CN were taken from Perez-Beaupuits et 
al. ( 20071 1. The peak main beam temperatures of the lower tran- 
sitions of HCO^ were provided by M. Krips (private communi- 
cation; Krips et al. 2008). 

The line ratios were computed considering uncertainties of 
10% in the reported beam efficiencies ti,„i,, 5% in the main beam 
0„,h, and 10% in the estimated source size 0s . As described in the 
next section, 9s was assumed to be 1 .5" and the same for all the 
lines. Note that the high uncertainties given in Table |4] should 
be considered lower limits, since comparison of data observed 
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Table 5. Physical conditions estimated from line ratios" and line intensities. 



Molecule/ratio 


n(H2) 


Tk 
[K] 


[cmT^ km"' s] 




HCN 


10^ - lO'' 


> 20 


1013.2 _ 1015.2 


10-3° < Ti_o < 10^^ 
lO^^ < T3_2 < 10"** 
ICP'- < T4-3 < 10^** 


HNC 


10'' - lO'' 


10-90 


1012.8 _ iQlS.O 


10-' " < Ti_o < 10^ " 

1 < T,_2 < 10^ ** 
10-" « < T4_3 < 10^4 


CN 


10^ - - 10' 


10-20 


1014 _ 10152 


10-20 < T,3,,_o,,, < 10^^ 

1 < T2„,-l„, < 10^* 
10-° ^ < T3,,,_2,„ < 10' ^ 


HCO+I^ 


10^ ' - 10' 


> 30 


10" " -10'" 


lO-*" < T3-2 < 10-' 2 
10-2 2 < T4-3 < 10-' 2 



a) Results from the single-phase model with a source size 6^ = 1.5". 

b) In order to get the actual column density and optical depth, the values reported here must be multiplied by a pertinent line width, i.e. thermal 
plus turbulent width (FWHM) of the lines. According to Schinnerer et al. (2000), the velocity dispersion at the nuclear region of NGC 1068 is 
about 23 km s-'. 



with different telescopes may introduce additional uncertainties 
due to the presence of different unknown systematics between 
the telescopes. 

4.3. Modeling the excitation conditions 

We explore a wide range of possible excitation conditions that 
can lead to the observed line intensities and line ratios, in a sim- 
ilar way as described in Perez-Beaupuits et al. (I2007I I. In the 
radiative transfer models, the kinetic temperature Tk and the 
hydrogen density n(H2) can be constrained from the line ratios 
(Table|4]l, whereas the column density (per line width) is mostly 
constrained by the line intensities. We used the intensity of the 
J=3-2 line of HCN, HNC, and HCO+, and the Nj = 25/2 - I3/2 
line of CN, since these are common in the ratios between the 
low-7 and high-7 lines. The intensity of these lines was cor- 
rected by beam dilution as described below. We use the radiative 
transfer code RADEX to build a cube of kinetic temperature Tk, 
number density ^(Ha) and column density per line width N/AV 
for each molecule. 

For most of the molecules we explored collisions only with 
H2, since this is considered to be the most important astrophysi- 
cal process. Other collision partners can be H and He. Although 
their collision cross sections are comparable, H2 is about 5 times 
more abundant than He, and H is at least one order of magnitude 
less abundant than H2 in the dense cores of molecular clouds 
(e.g., Meijerink & Spaans 2005). Hence, the effects of collisions 
with these other partners can be neglected. Nevertheless, in con- 
trast with HCN, HNC and HCO+, the open shell molecular struc- 
ture of the CN radical makes it more interactive with electrons. 
Therefore, we also explore the CN-g- collision as an additional 
process that could aid the excitation of CN, since the electron 
abundance can be enhanced by a high ionization degree in PDRs 
and XDRs. 

The volume density explored ranges between 10"* cm-"* and 
10' cm--^, the kinetic temperature varies from 4 K to 200 
K, and the column density per line width lies between 10'*' 
cm^^km"' s and 10'^ cm-^km-' s. In order to get the actual 
column density and optical depth, the values that will be re- 



ported here must be multiplied by a line width. According to 
Schinnerer et al. (2000), the velocity dispersion at the nuclear re- 
gion of NGC 1068 is about 23 km s"' . So, we could expect local 
velocity dispersions (line widths) ranging between 3 km s"' and 
lOkms-'. 

The original RADEX code was modified in order to include 
dust background emission as a diluted blackbody radiation field, 
in the same way as done by Poelman & Spaans (2005). We con- 
sider the emission component at 34 K described in Spinoglio 
et al. (2005), scaled by the continuum optical depth Te{v) at a 
particular frequency. Since in our models we consider up to ten 
transitions, we estimated rdv) as increasing linearly with the fre- 
quency of the transition (from Fig.4c in Spinoglio et al. 2005). 
We modeled the total background radiation as a composite be- 
tween the Cosmic Background Radiation (CMB), as a blackbody 
function at 2.73 K, and the diluted dust radiation, which was es- 
timated as Tc X B{Tdust), where B{Tdust) is the Planck function at 
Tdust = 34 K and rdv) - Tioo;im(100//m//l). From Hollenbach et 
al. (1991) we adopted the value rioo/jm = 10"^, hence the con- 
tinuum optical depth becomes Tc(v) - 10-^[cm]v/c, where v is 
the rest frequency of the transition in [Hz], and c is the speed of 
light in [cm/s]. 

The inclusion of the diluted continuum emission does not af- 
fect the results found for HCO^ but decreases (from 10** 'cm-^ 
to 10^ ■^cm--') the minimum density at which solutions are found 
for CN, and it increases (by ~ 10 - 20 K) the minimum tem- 
peratures found for HCN and HNC, at the lowest densities, in 
comparison with a model that considers only the CMB as back- 
ground radiation. The lower limits of the column densities per 
line width increases as well. 

The physical conditions of HCN, HCO+, and HNC were 
modeled using the collisional data available in the LAMDAQ 
database (Schoier et al. 2005). In the case of CN, we used the 
collision rates estimated by Fuente et al. (1995). We do not con- 
sider hyperfine splitting of CN, HNC and HCN, since they are 
not resolved in extragalactic systems. 

In Perez-Beaupuits et al. ( |2007l l the excitation conditions 
of HCN and HNC overlapped, assuming that their emission 

3 http://www.strw.leidenuniv.nl/~moldata/ 
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emerges from the same region. Since we now have more tran- 
sitions, we can constrain the temperature and density for each 
molecule separately, and then see if they agree or not. 

Once the starburst contribution to the low-7 lines is removed, 
the bulk of the emission for all the dense gas tracers mostly arises 
from the CND. There should not be a significant difference in the 
sizes of the CND as seen in HCN, HCO+, HNC and CN, if all the 
lines can be modeled by a single phase. So we assume the same 
source size for all the J lines to correct them for beam dilution. 
From the interferometer maps of CN, HCN and HCO^ available 
in the Uterature for NGC 1068 (e.g. Tacconi et al. 1994, Kohno et 
al. 2001, Garcia-Burillo et al. 2008), we estimate that the source 
size (FWHM) of these molecules ranges between about 1 and 2 
arcseconds. 

We can model the physical conditions for most of the line 
ratios of most of the molecules, considering an average source 
size 0s = 1.5". A summary of the estimated physical conditions 
is provided in Table|5] and a summary of the main results is pre- 
sented in the next section. Details of the modeling and analysis 
of the uncertainties are given in the Appendix. 

4.4. Physical environment and abundance ratios 

All the solutions for HCN are found for temperatures higher 
than 20 K, with a degeneracy between the kinetic temperature 
and molecular hydrogen density. The temperatures at which so- 
lutions can be found for HNC are limited up to about 90 K. 

The solutions found for both CN J^j and line ratios over- 
lap in an small region, with a narrow temperature range (between 
10 K and 20 K). However, if we use an smaller source size the 
temperature range can go from 20 K up to 100 K. The analy- 
sis of the uncertainties in CN is presented in Appendix B.l. We 
found that the effect of electrons as a secondary collision partner 
of CN, is not important in a PDR environment due to the rela- 
tively low electron abundance (~ 10 "*) with respect to H2 and 
because the abundance of CN is very low where n(e") is consid- 
erably higher (at the outer layers of a molecular cloud, close to 
the ionization front). In an XDR environment, small effects can 
be expected from collisions with electrons, mainly for the high-7 
lines. 

Using 0s - 1-5" our model is not able to reproduce the 
HCO^ line ratio and the HCO^ 7=3-2 line intensity simulta- 
neously. But it is able to reproduce the HCO^ However, our 
model can reproduce the observed HCO^ |^ ratio, and HCO^ 
7=3-2 intensity, if we assume either a larger source size, or a 
larger starburst contribution factor Detailed analysis of these al- 
ternatives is described in Appendix B.2. 

We estimate an A^(HNC)/A^(HCN) column density (abun- 
dance) ratio lower than unity in the CND of NGC 1068. This 
result cannot be interpreted as a signature of a pure PDR or XDR 
environment because the HNC abundance can be decreased 
due to high temperatures produced inside a molecular cloud by 
mechanisms like turbulence and shocks (Loenen et al. .2008) . 
The high SiO abundance observed in the CND of NGC 1068 can 
be a direct evidence of the possible contributions from mechani- 
cal heating (shocks) and dust grain chemistry induced by X-rays 
(Garcia-Burillo et al. 2008). Hence, a more elaborate PDR/XDR 
model that includes both mechanical heating and grain surface 
chemistry is needed for the interpretation of our results. 

The mean A^(CN)/A^(HCN) column density ratios (0.7-0.9) 
estimated with 6s = 1.5" can be easily found in an XDR 
environment, but the predominance of this component cannot 
be concluded from the A^(CN)/A^(HCN) ratio only, since ratios 



~ 1 .0 are also expected in a PDR component (Lepp & Dalgarno 
1996, Meijerink et al. 12007 ). The mean column density ratios 
2 < A^(CN)/A^(HCN) < 4 estimated with 63 = 1", are tenta- 
tively more consistent with an XDR/AGN environments (Lepp 
& Dalgarno 1996, Meijerink et al. 120071 ). However, these and 
higher [CN]/[HCN] abundance ratios have also been found in 
PDR/starburst scenarios (e.g. Fuente et al. 2005). The fact that 
we find a relatively low A^(CN)/A^(HCN) column density ratio, 
with respect to what would be expected in a pure XDR sce- 
nario, could be explained by an overabundance of HCN due to 
the grain-surface chemistry suggested by Garcia-Burillo et al. 
(2008). 

Our model indicates that the emission from the high-/ HCO^ 
lines emerge from gas that does not co-exist with HCN, HNC 
and CN, in the nuclear region of NGC 1068. The warmer and 
denser gas traced by the high-7 lines is only an small fraction 
(0.5% - 10%) of the total HCO+ gas. Most of it is confined to 
the lower transitions. Details of the analysis and interpretation 
of the abundance ratios are presented in Appendix C. 



5. Conclusions 

We have used HARP on JCMT to observe the 7=4-3 line of 
HCN, HNC, HCO+ , and the A^y = 25/2 - I3/2 and A^y = 35/2 - 
2g/2 transitions of CN, in the nuclear region of NGC 1068. We 
estimated the excitation conditions and abundance ratios of these 
molecules in the CND from radiative transfer models, assuming 
a single phase gas for all the molecules and lines. We compared 
column density ratios with predictions from PDR/starburst and 
XDR/AGN models. 

We deduced a first-order estimate of the starburst contribu- 
tion to the 7=1-0 lines, out of independent observations with 
different telescopes. The starburst contribution of the higher-7 
lines should be negligible because of the smaller beam sizes at 
the higher observing frequencies. 

We first expected, and assumed, that the source sizes, seen 
by the different molecules and transitions, should not be signifi- 
cantly different. We estimated the source sizes for HCN, CN and 
HCO^ from the different interferometer maps available in the 
literature, and used a common value of 1 .5" for all the observed 
lines. The results of our models showed that the HCN, HNC and 
partly CN, can co-exist in gas with similar kinetic temperature 
(10-20 K) and hydrogen density (lO**-^^"^-^ cm^^). 

However, the analysis of the uncertainties in source size, 
as well as in the starburst contribution factor, shows that the 
HCO^ molecule seems to trace a totally different gas phase, 
both with respect to the other molecules and its own different 
transitions. The HCO^ |5| line ratio indicates that the emis- 
sion of these lines emerges from gas with densities larger than 
10^^ cm"-' and temperatures higher than 30 K. Instead, the 
HCO^ line ratio would traces less dense (n(H2) < 10^ cm"^) 
and colder (Tk < 20 K) gas, when using models with either an 
smaller source size or a higher starburst contribution factor 

The physical conditions estimated from the line ratios are 
based on the assumption of a single-phase model with the same 
size of the emitting region. However, a single-phase scenario 
clearly does not hold for the HCO^ molecule. The lack of an 
independent measurement of the HCO^ 7=1-0 line, and the as- 
sumption of the same starburst contamination factor for this line 
as the one found for HCN from the IRAM 30m data, introduces 
an additional uncertainty. Nevertheless, the fact that the high- 
7 lines of HCO^ clearly trace higher densities and temperatures, 
with lower HCO^ column densities, indicates that this gas is less 
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extended throughout the CND, in terms of volume and maybe in 
surface as well. 

The Une intensity ratios of CN/HCN and HNC/HCN de- 
crease with increasing rotational quantum number J. The 
HCO^/HCN intensity ratios, instead, do not follow this trend. 
The HCO"^ /=4-3 line intensity, compared with the lower tran- 
sition lines and with the HCN 7=4-3 line, support the influence 
of a local XDR environment. 

TheA^(CN)/A^(HCN) ~ 1 -4 column density ratios estimated 
with the single-phase model are indicative of an XDR/AGN en- 
vironment with a possible contribution of grain-surface chem- 
istry induced by X-rays or shocks. 

The single-phase model, with a common source size seen by 
aU the molecules and transitions, works well for HCN and HNC. 
However, HCO^ shows clear evidence that different lines can 
trace completely difl'erent gas phases. Besides, assuming an av- 
erage source size for all the molecules and lines can be a good 
first approximation, but has the clear drawback that estimates de- 
rive mostly from the /=l-0 interferometer maps, which do not 
trace very warm and dense gas, and are sensitive to beam smear- 
ing efl'ects. High resolution maps of the higher-7 lines may show 
slight differences in source sizes that can be considered negligi- 
ble, at a first glance. But, as was shown with CN and HCO^, even 
an error of a fraction of an arcsecond in the estimated source size 
can make a big difference in the physical conditions derived from 
the radiative transfer models. High resolution maps of the difl'er- 
ent transitions used in the models would aid to reduce the uncer- 
tainties in the models. In addition, a more elaborated PDR/XDR 
model that includes mechanical heating and grain-surface chem- 
istry, will be necessary to properly account for the different con- 
tributing factors. 

The spatial extent and morphology of the emitting regions of 
different J lines, as well as the detailed dynamics and physical 
environment of the high density tracers, can be tested with the 
higher spatial resolution maps that ALMA wiU provide in the 
future. 
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Appendix A: Modeling the physical conditions of 
high density tracers 

A.1. The physical environment of HCN 



The top panel of Figure lA.ll shows the mean column density 
per hne width A^(HCN)/Ay cm"^ km ' s (the mean between the 
maximum and minimum N/Av that yield valid solutions), for all 
the explored densities n(H2) and temperatures Tk that can repro- 
duce, within lo", the observed HCN 2^ and HCN |^ line ratios 
and the intensity of the HCN 7=3-2 line. The bottom panel of 
Figure [ATI shows how A^(HCN)/Al' cm^^km"' s changes with 
n(H2) at different temperatures. These curves are easier to com- 
pare with the output of Large Velocity Gradient (LVG) models 
(e.g. Goldreich & Scoville 1976) commonly found in the litera- 
ture. To make these curves more clear we do not show the error 
bars (or uncertainties) of A^(HCN)/Al' for each temperature in 
the plot. But they range between 25.0% - 48.8% at 30 K, 24.2% 
- 40.6% at 50 K, 18.4% - 38.0% at 70 K, and 21.2% - 36.4% 
at 90 K. Note that the uncertainties decrease towards the higher 
n(H2) densities. 

The lower and upper limits of the HCN column density per 

^13.2 ^^-2 g ^j. 



line width A^(HCN)/Al' cm "^km s are 10 cm "^km 



^ km"' I 

the highest density of about 10^ cm"-', and lO'^ "* cm"^ km"' s at 
the lowest density explored of 10** cm"^, respectively. There is 
also a narrower temperature region (around 20 K) with densities 
higher than 10* cm"^, where solutions for the observed ratios 
and intensities are also possible. The optical depths r of each 
line are summarized in Table |5] 

All the solutions are found for temperatures higher than 20 
K, with a clear degeneracy between the kinetic temperature and 
molecular hydrogen density. That is, a given column density can 
be obtained with either high and low «(H2), or low Tk and 
high n(H2). 

A.2. The physical environment of HNC 

In the case of the double peak structure observed in the lower 
transition lines of HNC, only the main peak (at velocity ~ 
1073 ±13 kms"') was considered in the analysis, since this is 
the component that is closer to the central velocity observed in 
the HNC 7=4-3 line. 

The top panel of Figure IA.2I shows all the excitation con- 
ditions for which the observed HNC |5| and HNC |5| Une 
ratios and the intensity of the HNC 7=3-2 line can be repro- 
duced within Icr. The mean column density of HNC per line 
width (NiHNC)/ Av cm"^ km"' s) is shown in the contour plot. 
The bottom panel of Figure lA. 21 shows how the column density 
changes with n(H2) at different temperatures. The corresponding 
eiTor bars ranges ai-e 6.4% - 35.9% at 30 K, 7.6% - 32.3% at 50 
K, and 8.3% - 26.0% at 70 K. The error bars decrease towards 
the higher n(H2) densities. 

The estimated A^(HNC)/Al' cm"^km"'s column den- 
sity per line width ranges between 10'^'^ cm"^km"'s and 
lO'^"' cm"^km"' s. The n(H2) densities required to reproduce 
these ratios range from 10"* cm"^ (and probably lower than that, 
when considering higher temperatures outside of our explored 
grid) and 10^ cm""* at temperatures of about 10 K. A summary 
of the corresponding optical depth t can be found in Table |5] 

In contrast with HCN, the temperatures at which solutions 
can be found for HNC are limited up to about 90 K, for the 
lowest densities explored. At a density of 10^^ cm"^, Tk ranges 
between 10 K and 70 K. But at densities > 10* cm"-' only tem- 
peratures lower than 30 K are allowed. 
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Fig.A.l. Top - Excitation conditions modeled for 
and 1^ line ratios of HCN. The contour lines corre- 
spond to the mean column density of HCN per line width 
(A^(HCN)/Ay cm"2km"' s), in the region where the estimated 
line ratios and the 7=3-2 line intensity reproduce, within Icr, 
the observed values. Bottom - Transversal cuts of the column 
density per line width at different temperatures modeled above. 
The column densities are of the order of 10'"' cm"^km"' s, in 
linear scale. 



A.3. The physical environment of HCO^ 

Using the common source size of 1.5" and the starburst con- 
tribution factor of 0.45 in the HCO^ 7=1-0 line, our model 
is not able to reproduce the observed HCO^ |5| line ratio 
and HCO^ 7=3-2 line intensity. The top panel of Figure IA.3I 
shows the mean column density of HCO^per line width 
(A^(HCO+)/Ay cm"2km"' s) modeled for the HCO+ |^ Hne 
ratio and HCO^ 7=3-2 line intensity. The bottom panel of 
Figure |A3] shows how the column density changes with n(H2) at 
different temperatures. The corresponding error bars ranges are 
29.9% - 30.1% at 50 K, and 29.3% - 30.0% at 70 K, and 29.0% 
- 30. 1 % at 90 K. In this case the error bars decrease towards the 
lower n(H2) densities. 

The model shows that the high-7 lines trace gas with den- 
sities larger than 10^^ cm"^^ and temperatures larger than 30 
K. Solutions with temperatures lower than 30 K could also 
be found at densities larger than those explored in this work 
(n(H2) > 10^ cm"^). The columns range from lO"-' cm"^ km"' s 
to 10'^^ cm"^km"'s, and both lines are optically thin (t < 
10"'^) over the whole range of columns (Table |5]l. The analy- 
sis of the uncertainties in the HCO^ model is discussed in the 
Appendix B.2. 
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Fig.A.2. Top - Excitation conditions modeled for the 
and line ratios of HNC. The contour lines corre- 

spond to the mean column density of HNC per line width 
(N(HNC)/Av cm-^km"' s), in the region where the estimated 
line ratios and the 7=3-2 line intensity reproduce, within Icr, 
the observed values. Bottom - Transversal cuts of the column 
density per line width at different temperatures modeled above. 
The column densities are of the order of 10'"* cm -km"' s, in 
linear scale. 



A.4. The physical environment of CN 

Only the main spingroups of each transition were considered for 
modelling the physical conditions of CN. The hyperfine struc- 
ture of CN is not included in the model. Hence, the radiative 
lines of CN are described just by the quantum numbers and J. 
The collisional data are the same as used by Fuente et al. (I1995I I. 

The excitation conditions derived from the two line ratios 
with the one-phase model (with a source size of 1.5") are shown 
in the top panel of Figure IA.4I In contrast with HCN and HNC, 
these conditions overlap just in an small region, with a narrow 
temperature range. At the lowest density for which solutions are 
found (n(H2) = 10^^ cm"^), the kinetic temperature of the gas 
Tk is between 14 K and 16 K. At higher densities (n(H2) > 
10^ cm"^) the gas can be as cold as 10 K. 

The bottom panel of Figure IA.4I shows how the col- 
umn density changes with n(H2) at Tk = 14 K. The 
dashed lines correspond to the upper and lower limits of 
A?(HCO+)/Ai/ cm^^km"' s. Note that for n(H2) = 10^^ cm"^ 
and n(H2) = 10^ cm"^ the limits converge, which means there is 
only one solution at T^: = 14 K for those densities. 

The column densities per line width have a lower limit of 
about 10'** cm"- km"' s at n(H2) = 10^ cm"-' and an upper limit 
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Fig.A.3. Top - Excitation conditions modeled for the |^ line 
ratio of HCO^ . The contour lines correspond to the mean column 
density of HCO+per Hne width (A^(HCO+)/A(7 cm"- km"' s), in 
the region where the estimated line ratios and the J-3-2 line 
intensity reproduce, within Icr, the observed values. Bottom - 
Transversal cuts at different temperatures of the column density 
per line width modeled above. The column densities are of the 
order of 10" cm"- km"' s, in linear scale. 

of 10'^^ cm"^ km"' s at a density of 10'*^ cm"^. The correspond- 
ing optical depths of each line can be found in Table |5] 

Besides molecular hydrogen, we also explore the effects of 
electrons as a second collision partner for CN. These could be 
significant both in PDR and XDR environments, due to enhance- 
ment of the ionization degree by radiation. The collision rates are 
the same as used in Black & Van Dishoeck ( 119911 1, and were ob- 
tained from J.H. Black (private communication). These rates are 
available only for Tk = 20 K. 

In a PDR environment, the CN emissivity peaks at a total 
hydrogen column density of about 10^'^ - 10-^ cm"^. At these 
depths the electron abundance is ~ 10"^, and the CN column 
density is of the order of 10'"* cm"^. We found that for densi- 
ties n(H2) > 10^ cm""* the effect of collisions between CN and 
electrons is negligible. 

The only region where the electron abundance can be larger 
is at the edge of a PDR cloud, that is A^h ^ 10^'^ cm"^. There 
the electron abundance is still about four orders of magnitude 
lower than the total hydrogen density, but it can be at least two 
orders of magnitude larger than the H2 density (Meijerink & 
Spaans 120051 ). At those shallow depths, however, the column of 
CN is not significant (A^(CN) < lO"* cm"^). In order to boost the 
electron abundance to higher levels, in the region where most of 
the CN emission originates, we would require a gas phase car- 
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Fig.A.4. Top - Excitation conditions modeled for the J'' 
and J'^ line ratios of CN. The contour lines corre- 

•^5/2-13/2 

spond to the mean column density of CN per line width 
(N(CN) / Av cm"^ km"' s) in the region where the estimated line 
ratios and the Nj - 2^/2 - I3/2 line intensity reproduce, within 
Icr, the observed values. Bottom - Transversal cut at Tk = 14 
K of the column density per line width modeled above. The 
values are of the order of 10'^ cm^^km"' s, in linear scale. 
The dashed lines correspond to the upper and lower limits of 
A^(CN)/Aiy cm-^km"' s. 

bon abundance of about 3-4 times Solar. However, these higher 
abundances are not supported by other works (e.g. Kraemer et 
al.i998) 

In an XDR, the ambient conditions along the cloud are differ- 
ent than those found in a PDR, with larger (10^^ - 10 "*) relative 
electron abundances (Meijerink & Spaans I2005I) . In the main 
emitting region (A^(H) > 10^^ cm"^) the electron density n(e^) 
is expected to be about 10 cm"-' if n(H2) ~ 10^ cm"-'. These 
densities produce changes of the order of 10% in the excitation 
temperatures ^ of CN, with respect to those obtained when us- 
ing only H2 as collision partner, and the column density N(CN) 
has to decrease with about 50% to get the same line strengths. 

Hence, the effect of electrons as secondary collision partner 
of CN, is not important in a PDR environment. In an XDR envi- 
ronment, small effects can be expected. 

Appendix B: Analysis of the uncertainties 

We analyze here how the uncertainties in the input parameters 
of our models (4-3/3-2, 3-2/1-0 line ratios and the J-3-2 line 
intensities) propagates to the solutions we find. The two main 
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Fig.B.l. Excitation conditions modeled for the l?'' and 

3 —2 

J^'' line ratios of CN, using a source size of 1". The con- 

^5/2^13/2 

tour lines are as defined before. Considering this smaller source 
size, the range of temperatures where solutions are found can go 
up to 100 K, and densities can be as low as 10"*-^ cm"-'. 



uncertainties in our models are the estimated common source 
size (Os - 1.5") for all the molecules and transitions, and the 
first order estimate of the starburst contribution factor /sb for 
the 7=1-0 lines. The starburst contribution factor affects only 
the line ratios which in turn will modify, to some extent, 
the combination of temperatures (Tk), densities (n(H2)) and col- 
umn densities (A^) for which solutions are found. On the other 
hand, the source size affects mostly the 7=3-2 line intensities 
(used to constrain the radiative transfer models) which in turn 
affect mostly the range of temperature and column densities of 
the solutions. The source size is also present in the 4-3/3-2 and 
3-2/1-0 line ratios, as we correct these lines for beam dilution. 
However, because the source size is about one order of magni- 
tude smaller than the size of the respective beams, its effect in 
the line ratios is negligible. 

The solutions for the HCN and HNC line ratios are less sen- 
sitive to these two parameters since the range of temperatures 
and densities for which solutions exist is large enough (which 
allows for more flexibility in the Tk vs n(H2) space) and be- 
cause the 1^ line ratios, which basically define and constrain the 
overlap with the solutions found for the ratio between the lower- 
7 lines, are independent of /sb . The effect is reflected mostly in 
the column densities per line width due to changes in the source 
size. 

B. 1 . Uncertainties in CN 

The solutions for CN are particularly sensitive to the source size. 
As mentioned before, the starburst contribution factor affects 
only the ratio between the lower 7-lines. The range of Tk and 

n(H2) for which we find solutions for the CN l"' line ratio 

25/2-I3/2 

is smaller than that for the ratio. Hence, the final solu- 

13/2-U1/2 

tions for CN (given by the overlap between the solutions found 
separately for the ratios between the low-7 lines and the high-7 

lines) is constrained by the solutions found for the J^'' line 

25/2-13/2 

ratio. 

Because the solutions found for the CN are already con- 
strained, little changes in the source size have a larger impact 
than for HCN and HNC. If we assume a larger source size of 2", 
the solutions for CN are restricted to an small range of temper- 
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ature around 10 + 4 K and for densities between 10^^ cm"^ and 
10^^ cm On the other hand, if we assume an smaller source 
size of 1", the solutions for CN line ratios can be found in a 
larger region of Tk and n(H2) than those found for a source size 
of 1.5", as shown in Figure IBTTI The smaller source size of 1" 
increases the estimated CN A^y = 2g/2 - I3/2 radiation temper- 
ature by a factor 2.25, which is reflected mostly in the range 
of kinetic temperatures at which we can find solutions for the 
observed line ratios. When using 9s = 1.5" the maximum tem- 
perature where solutions can be found is 20 K, at a density of 



10 cm" (Figure |A~4| |. But when using an smaller source size 
(Figure IbTI i the temperature range can go from 20 K up to 100 
K, at the same density. 



B.2. Uncertainties in HCO+ 

On the other hand, the HCO+ |5| is aff'ected by both, 0s and 
/sB, and the fact that we do not find solutions for this ratio may 
be due to the uncertainties in these two parameters. We explored 
different alternatives and we found that our model can reproduce 
the observed ratio, and HCO^ J -3-2 line intensity, if we assume 
either a larger source size {6s - 2"), or a larger starburst contri- 
bution factor (/sB = 0.60). The results are shown in the top and 
middle panels of Figure IBT2] A larger starburst contribution fac- 
tor increases the HCO^ |^ line ratio, and allows for solutions at 
slightly higher temperatures than assuming a larger source size. 
Note that in both cases solutions for densities < 10"* cm"^ are 
also possible. This result, however, rises a new question since 
these solutions do not overlap with the solutions found for the 
HCO+ 1^ ratio (Figure|A31l- 

The fact that the ratio between the higher transition lines 
traces denser (n(H2) > 10^'^ cm"-') and warmer (Tk > 30 K) 
gas than what the lower transition lines would indicate, and that 
they do not overlap either with the solutions found for HCN, 
HNC or CN, may raise some skepticism about the HCO^ 7=4-3 
line. So we also explored the results of our model, assuming that 
the HCO^ 7=4-3 line could be somehow overestimated (e.g., 
undetected calibration problems, or intrinsic instrumental differ- 
ences, like sensitivities, between the JCMT and the IRAM 30m 
telescopes). 

In section 2 we mentioned that HARP observations of the 
calibration source showed less flux than the reference spectra 
observed with the former receiver B3. Those differences account 
for factors between 1.1 and 1.5, that where used to correct each 
of our scans before adding them up. So, if we reduce the final 
HCO"^ 7=4-3 line by a factor 1.5 (which would be the worst 
case scenario) the HCO^ |^ line ratio reduces to 0.91 + 0.31 
and the new solutions would be as shown in the bottom panel of 
Figure IBT2] It can be seen that, even in this worst case scenario, 
there would not be an overlap between the solutions found for 
the ratio between the lower 7 lines and those found for the higher 
7 lines. 

This result would imply that the lower and higher 7-lines of 
HCO^ trace different gas phases, and a single-phase model may 
not be the most appropriate to reproduce the observed HCO^ ra- 
tios and intensities. Hence, a two-phase model, where different 
sizes of the emitting region could be seen by different 7 lines, 
may be a better approach for HCO^. In fact, Krips et al. (2008) 
mentioned that the size of the emitting region decreases with 
increasing 7 line, according to recent SMA data. A multiple- 
phase model can also be applied to the HCN, HNC and CN 
molecules. Among molecules, those which are easily dissoci- 
ated should have small sizes. Hence, a decreasing source size 
with increasing 7 line seems natural. However, the only way to 
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Fig. B.2. Top - Excitation conditions modeled for the jf^ line 
ratio of HCO^, when considering the estimated starburst con- 
tribution factor fsB - 0.45, but a larger source size 9s - 2". 
Middle - Excitation conditions modeled for the line ratio of 
HCO^, when considering a larger starburst contribution factor 
fsB - 0.60, and the common source size 9s - 1.5". Bottom - 
Excitation conditions modeled for the |5| lin^ ratio of HCO^, 
when reducing the HCO^ 7=4-3 line intensity by a factor 1.5. 
The higher transition 7=4-3 still traces mostly the dense and 
warm gas (n(H2) > lO^"* cm""*, T^ > 20 K), whereas the 
1^ line ratio indicates that those lines would trace less dense 
(n(H2) < 10^-2 cm-^) and cold (10 K < < 30 K) gas. 



constrain these models would be by using interferometer maps 
(fluxes, beam deconvolution, etc.) to accurately estimate the ac- 
tual source size seen by the different 7 lines. Nevertheless, this 
is a task that is beyond the scope of this work. 
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Table C.l. Line intensity ratios between molecules. 



Molecules 






Transition" 






y= 


1-0 




3-2 


y=4-3 


HNC 
HCN 


0.71 


±0.45 


0.26 


±0.08 


0.16+0.05 


CN 
HCN 


0.89 


±0.69 


0.54 


±0.15 


0.14±0.04 


HCO+ 
HCN 


1.10 


±0.59 


0.13 


±0.04 


0.36±0.10 



a) The lines correspond to J=l~0, 7=3-2 and 7=4-3, for HCN, 
HNC and HCO^. In the case of CN, the lines are actually Nj = 
I3/2 - O1/2, Nj = 2s/2 - I3/2 and Nj = 35/2 - 25/2, respectively. 



Appendix C: Line intensity and abundance ratios 

The line intensity ratios, with respect to HCN, are summarized 
in Table IC.ll The intensities used correspond to the peak an- 
tenna temperature of the main component of the gaussian fits, 
corrected for starburst contribution, beam efficiency, and beam 
dilution, assuming a source size of 1.5" for all the lines. We also 
find that most of the line intensity ratios (the HCO+/HCN/=4- 
3 ratio is the exception that is discussed below) decrease with 
increasing rotational quantum number 7, similar to what was 
found by Krips et al. (2008) for the HCN/CO ratio. This was 
also noticed for the HNC/HCN 7=1-0 and 7=3-2 line ratios 
observed in NGC 1068, and other Seyfert galaxies, by Perez- 
Beaupuits et al. (2007). Since the beam sizes are comparable, 
this could be an indication that the higher-7 levels of CN and 
HNC are less populated than those of HCN. 

This is what we would expect assuming collisional excitation 
of the molecules, where Tg^ is proportional to tx«(H2). From the 
excitation conditions modeled in Section §4.3, the HNC column 
densities, and hence thnc, tend to be lower than that of HCN. On 
the other hand, even if the column densities that we find for CN 
are as high (or higher in some cases) than the columns found for 
HCN, the optical depth tcn tends to be lower than thcn because 
of the higher (fine and hyperfine) splitting of the rotational levels 
of CN. Thus, the high energy levels of CN and HNC would be 
less populated than those of HCN. 

In the following sections we address the relative abundance 
issue by estimating the abundance ratio as the ratio between the 
column densities estimated from the radiative transfer models 
described before. 



C.1. HNC/HCN 

The range of temperatures where we can analyze the 
HNC/HCN ratio is limited by the solutions found for HNC, 
which go up to ~ 90 K at the lowest density explored (n(H2) = 
10^ cm""*). The top panel in Figure IC.ll shows the mean 
A^(HNC)/A^(HCN) column density ratio for the temperatures and 
densities where the solutions found for each molecule (top pan- 
els of Figs. IA.2I & IA.II 1 overlap. The mean column density ra- 
tio ranges between 0.10 and 0.18, with errors that vary between 
30% and 55% of the mean value. Note that similar ratios can be 
found at two different temperatures for a particular density. The 
bottom panel in Figure IC.ll shows the mean value, and corre- 
sponding upper and lower limits, of the A^(HNC)/A^(HCN) ratio 
at 20 K and 50 K. At these temperatures the ratios are quite sim- 
ilar (within 15%) for the density range n(H2) = lO^'^*"^'^ cm"-'. 

The fact that the line intensity ratios (Table IC.lb are also 
lower than unity, indicates that the bulk of the HNC and HCN 
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Fig. C.l. Top - Ratios between the mean column densities of 
HNC and HCN (A?(HNC)/N(HCN)). Bottom - Transversal cut 
of the ratios at different temperatures. The dashed and dashed- 
dot lines correspond to the upper and lower limits of the ratios at 
20 K and 50 K, respectively. 



emission emerges from warm gas (T^ > 30 K). This agrees with 
observations in the vicinity of the hot core of Orion KL, and ex- 
perimental and theoretical data, where the HNC/HCN line ratio 
decreases as the temperature and density increase (e.g. SchiUce 
et al. 1992; Talbi et al. 1996; Tachikawa et al. 2003). 

On the other hand, the ratios A^(HNC)/A^(HCN) < 1 esti- 
mated with our models cannot be directly interpreted as a sig- 
nature of a pure PDR or XDR environment in the CND of 
NGC 1068. The HNC abundance can be decreased due to tem- 
peratures higher than traditionally expected, produced deep in- 
side a molecular cloud by mechanisms other than radiation, like 
turbulence and shocks (Loenen et al. I2008I I. If the tempera- 
ture is higher than 100 K, the conversion of HNC into HCN 
is more efficient and HNC is suppressed (Schilke et al. 119921 
Talbi et al. 1 19961 1. These high temperatures are not found in re- 
gions where the abundance of HCN and HNC is high enough 
to be detected, for traditional PDR or XDR models (Meijerink 
& Spaans 2005). The high SiO abundance observed in the CND 
of NGC 1068 can be a direct evidence of the possible contri- 
butions from mechanical heating (shocks) and dust grain chem- 
istry induced by X-rays (Garcia-Burillo et al. 2008). Hence, a 
more elaborated PDR/XDR model that includes both mechani- 
cal heating and grain surface chemistry will be needed to further 
understand the results of our radiative transfer models. 
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Fig.C.2. Top - The A^(CN)/A^(HCN)) column density ratios for 
the small range of density at = 18 K, where the physical 
conditions found for CN and HCN overlap when using the one- 
phase model with 8^ - 1.5". Middle - Ratios between the mean 
column densities of CN and HCN, when using a source size Qs — 
1" for CN. Bottom - Transversal cuts, at different temperatures, 
of the ratios showed in the middle panel. The dotted, dashed- 
dotted and dashed lines correspond to the upper and lower limits 
of the ratios at 18 K, 30 K and 50 K, respectively. 



C.2. CN/HCN 

When using the one-phase model with a source size 6s - 
1.5", the physical conditions estimated for the CN and HCN 
molecules (top panel in Figs. |A.4l & |A.U overlap only in an small 
region around n(H2) ~ 10^'^ cm""* and Tk ~ 20 K. The top panel 
in Figure |C2] shows the mean A^(CN)/A^(HCN) column density 
ratio and the corresponding upper and lower limits. The maxi- 
mum ratio of 0.93 + 0.07 is found at n(H2) ~ 10^ '* cm"-'. The ra- 
tio decreases almost linearly with density (while the uncertainty 



increases), and it reaches the minimum ratio of 0.73 + 0.2, with a 
larger uncertainty of ~ 29%, at a density of n(H2) ~ 10^^ cm"-'. 
The ratio then increases slowly, with a constant uncertainty, up 
to 0.76 + 0.2 at n(H2) ~ 10^-^ cm^^ 

If we assume a larger source size 6s - 2" for CN, the 
solutions obtained do not overlap with those found for HCN. 
However, when using an smaller source size 6s - 1" for CN, 
there is a large overlap in the physical conditions found for these 
molecules (top panel in Figs. IB.ll & IA.ll i. The middle panel in 
Figure lC2] shows the mean A^(CN)/A^(HCN) column density ra- 
tio for the new model. The new mean column density ratios now 
range between ~ 2 and ~ 4, with errors that vary between 20% 
and 50% of the mean value. The maximum mean ratios (> 3.6) 
are found at a density of ~ 10^ cm"^ and temperatures between 
40 K and 60 K, while the minimum mean ratios (< 2.3) can be 
found at a larger density range of n(H2) ~ 10^"^ cm"-' but at 
lower temperatures (Tk < 30 K). 

The bottom panel in Figure |C2] shows the mean value, and 
corresponding upper and lower limits, of the A^(CN)/A^(HCN) 
ratio at 18 K, 30 K and 50 K. At = 18 K the new ratios (with 
uncertainties ranging from 23.6% to 26.5%) are larger (rang- 
ing between 3 and 4) and are found at a higher density range 
(«(H2) > lO*" *" cm"^) than in the 1.5" source size model. On the 
other hand, at densities between 10^ cm"^ and 10^ cm"^ the ra- 
tios found for 30 K and 50 K are very similar (within 10%), but 
the uncertainties of the ratios at T/f = 30 K are larger (63.6%- 
72.4%) than those for = 50 K (39.0%-43.7%). These larger 
uncertainties for the ratios at = 30 K imply that the cor- 
responding upper limits increase to values > 5 at a density 
«(H2) ~ 10^-'' cm""*. Note that the uncertainties at Tk - 30/r, 
and hence the upper limits, decrease for densities < lO''-' cm"^. 



The mean column density ratios estimated with 6s - 
1.5" can be easily found in an XDR environment, but the 
predominance of this component cannot be concluded from 
the A^(CN)/A^(HCN) ratio only, since ratios ~ 1.0 are also 
expected in a PDR component (Lepp & Dalgarno 1996, 
Meijerink et al. I2007I I. The mean column density ratios 2 < 
A^(CN)/A^(HCN) < 4 estimated with 6s = 1", are tenta- 
tively more consistent with an XDR/AGN environments (Lepp 
& Dalgarno 1996, Meijerink et al. |2007| l. However, these and 
higher [CN]/[HCN] abundance ratios have also been found in 
PDR/starburst scenarios (e.g. Fuente et al. 2005). The fact that 
we find a relatively low A^(CN)/A^(HCN) column density ratio, 
with respect to what would be expected in a pure XDR sce- 
nario, could be explained by an overabundance of HCN due to 
the grain-surface chemistry suggested by Garcia-Burillo et al. 
(2008). 

Nevertheless, assuming an smaller emitting region for CN 
than for HCN introduces a new question regarding the chem- 
istry/physics driving the formation (or destruction) of these two 
molecules. If this is the case, then we would need to explain why 
CN is absent in the hypothetically more extended region covered 
by HCN. Perhaps this scenario could also be explained by the 
possible contributions from mechanical heating and dust grain 
chemistry suggested in Loenen et al. (2008) and Garcia-Burillo 
et al. (2008). However, exploring this alternative would require 
high resolution maps of at least SiO, HCN, CN and HNC, in ad- 
dition to a composite mechanical heating. X-rays and dust-grain 
chemistry model, to properly account for the different contribut- 
ing scenarios. This is an study that can be addressed in a follow 
up work. 
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C.3. HCO+/HCN 

Figure IA.3I indicates that the emission from the high-7 
HCO^ Hnes emerge from gas that does not co-exist with 
HCN, HNC and CN, in the nuclear region of NGC 1068. The 
column density (10" '' - 10'^'^ cm^^km^'s) estimated from 
the HCO^I^I ratio, and the possible solutions found for the 
HCO'*' y5§ ratio (top panels of Figure IB.2l i. also indicates that 
the warmer and denser gas traced by the high-/ lines is only an 
small fraction (0.5% - 10%) of the total HCO+ gas. Most of it is 
confined to the lower transitions. 

The main reason for the lack of co-existence is the HCO^ 
7=4-3 line. The possible solutions found for the HCO^ ratio, 
considering the uncertainties, would be consistent, in terms of 
density, with the HCN and HNC molecules, albeit at somewhat 
lower temperature. However, if we consider a larger source size 
of about 2" (as shown in Section 4.3), the HCO^ 7=3-2 line 
intensity will decrease, and solutions for temperatures up to 30 
K (at densities of a few times lO'* cm"-') will be possible, and 
the solutions for the high-7 ratio will require just slightly lower 
(A^(HCO^) ~ io"-^"i2.o cm^^km"' s) column densities per Une 
width. Hence, the HCO^ 7=4-3 line seems to indicate a different 
gas phase. 

Krips et al. (|2008) found that the HCN/CO line intensity 
ratios decrease with increasing rotational quantum number 7, 
for AGN dominated galaxies, including NGC 1068. We find the 
same trend in the HNC/HCN and CN/HCN line intensity ratios 
(Table |CB- The HCO+/HCN ratio, however, defies this trend. 
Interestingly, the higher-7 levels of HCO^ may be more popu- 
lated than its lower levels due to a local X-ray source (Meijerink 
et al. 2007). This result is consistent with an XDR, given that 
in strongly irradiated dense XDRs the HCO^ column builds up 
with depth to high values before the HCN does. Hence, the col- 
umn weighted temperature of the HCO^ molecule is higher, 
from which the HCO+ 7=4-3 line benefits. The HCN behav- 
ior with depth is more gradual, avoiding the strong separation 
between the low and high-7 lines (e.g., figure 9 of Meijerink & 
Spaans 2005). 

The 7=4-3 line ratio between the peak intensities of HCN 
and HCO^ (the inverse value is shown in Table 6) is about ~ 2.7, 
and is consistent with the ratio reported by Kohno et al. (2001). 
Instead, the velocity-integrated intensity ratio ^^qqI^ 7 = 4-3 is 
~ 3.7 (from Table 3), which is larger than the ratio between the 
peak intensities due to the smaller line width of the HCO^ 7=4- 
3 line. This number is right above the maximum value shown in 
Fig.3 of Kohno (2005). Interestingly, this places NGC 1068 in 
a distinguished position within the Kohno diagram, among their 
pure AGNs. 



